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Introduction 



43 

f} It has been recognized that the turbulent cross hclicity (correlation between the velocity 

and magnetic-field fluctuations) can play an important role in several magnetohydrodynamic 

S(MHD) plasma phenomena such as the global magnetic-field generation, turbulence suppres- 
sion, etc. Despite its relevance to the cross-helicity evolution, little attention has been paid to 
the dissipation rate of the turbulent cross helicity, eyy I n this paper, we consider the model 
expression for the dissipation rate of the turbulent cross helicity. In addition to the algebraic 
model, an evolution equation of £y/ is proposed on the basis of the statistical analytical theory 
of inhomogeneous turbulence. A turbulence model with the modeling of ew is applied to the 
solar- wind turbulence. Numerical results on the large-scale evolution of the cross helicity is 
O compared with the satellite observations. It is shown that, as far as the solar-wind application 

is concerned, the simplest possible algebraic model for ew is sufficient for elucidating the 
large-scale spatial evolution of the solar-wind turbulence. Dependence of the cross-helicity 
_l evolution on the large-scale velocity structures such as velocity shear and flow expansion is 

*^? also discussed. 

1 ' Keywords: Magnetohydrodynamic turbulence; turbulence model; cross helicity; dissipation 

rate; solar wind; 

in 
> 

\o 
I> 

In the magnetohydrodynamic (MHD) turbulent flow at high magnetic Reynolds 
number {Rm 3> 1), magnetic fields are considered to be frozen in plasmas, and 
move with the flow, [lj In such a flow, the induced magnetic field is often much larger 
J> than the originally imposed field. Besides, MHD waves such as the Alfven wave 

are considered to exist ubiquitously. The cross helicity, defined by the correlation 
between the velocity u and magnetic field b, is a possible describer of such MHD 
turbulence properties. Actually, the magnetic-field generation due to the turbulent 
cross helicity has been investigated. (2HH] 

As is well known, the total amount of cross helicity J v u • hdV, as well as that of 
the MHD energy J v (u 2 + b 2 )/2 dV, is an inviscid invariant of the MHD equations. 
Because of this conservative property, the turbulent densities of the MHD energy 
and cross-helicity, K = (u /2 + b /2 )/2 and W = (u' • b'), may serve themselves as a 
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good measure for characterizing the statistical properties of MHD turbulence (u': 
velocity fluctuation, b': magnetic-field fluctuation, (•••): ensemble average). 

The evolution equations of K and W are similar in form, and their mathemat- 
ical structures are quite simple. The evolution of K and W are determined by 
three constitutes: the production, dissipation, and transport rates. Firstly, the pro- 
duction rate is expressed by the correlation of turbulence fields coupled with the 
mean-field inhomogeneity. We note that the production rate of turbulence quan- 
tities such as (u ; ■ b') can be expressed exactly in the same as the counterpart of 
mean- field quantities such as U-B, but with the opposite sign (U: the mean veloc- 
ity, B: the mean magnetic field). This makes our interpretation possible that the 
drain of mean-field quantity gives rise to the generation of turbulence counterpart. 
So, the production rate represents how a quantity is supplied to turbulence by 
way of its cascading process (See Appendix |A|). Secondly, the dissipation rates of 
the turbulent MHD energy and cross helicity, e and ew, whose definitions will be 
given shortly in Section [2[ represent the effects of molecular viscosity and magnetic 
diffusivity coupled with the small-scale fluctuations. However, we stress the follow- 
ing point. The dissipation rates of the turbulent MHD energy and cross helicity, e 
and ew, can be considered from another aspect. In the intermediate range of tur- 
bulence, called the inertial range, the energy and cross helicity supplied from the 
energy-containing range compensate the energy and cross-helicity lost in the dissi- 
pation range. For this cascade picture of turbulence, the energy and cross-helicity 
transfer from lower to higher wavenumber ranges are most important quantities. 
In an equilibrium turbulence, e and ew represent these transfer rates of K and 
W, respectively. This makes the construction of the e and ew equation possible as 
we show later. Finally, the transport rates express the flux of a quantity that en- 
ters the fluid volume through the boundary. The expression for the transport rates 
suggests in what situation the quantities considered can be supplied to turbulence. 

Thanks to these clear-cut pictures associated with the evolution equation, the 
cross helicity (density) W, as well as the turbulent MHD energy (density) K, 
may play an important role in the turbulence modeling of MHD fluids. However, 
as compared with K and other pseudoscalar turbulence quantities such as the 
turbulence kinetic and magnetic helicities, only a limited attention so far has been 
paid to the cross helicity. 

In the context of homogeneous isotropic MHD turbulence, some important in- 
vestigations have been made on the decaying rate of the cross helicity or ew- It 
was shown that if there is a prevailed sign of the cross helicity in the initial state, 
the system goes towards a dynamically aligned state. The cross helicity scaled by 
the MHD energy grows towards +1 or -1 depending on the initially prevailed sign 
of the cross helicity. [T0HT3] 

In the context of inhomogeneous MHD turbulence, the cross helicity has been 
investigated mostly in the solar-wind research. By using spacecraft observations, 
detailed spectra of cross helicity have been examined. |14ffT6] In order to explain 
the large-scale behavior of the solar- wind turbulence, several models have been 
proposed. [TBI [T7j However, investigations related to the cross helicity are mostly 
concentrated on arguments of its production rate, and effects of large-scale inho- 
mogeneities such as the mean velocity shear have been discussed. Matthaeus and 
coworkers have employed a kind of algebraic model for the cross-helicity dissipa- 
tion rate. [T7ti20] Adopting this algebraic model of ew-, Usmanov et al. have recently 
performed a series of elaborated numerical simulations on the large-scale evolution 
of solar- wind turbulence. [21] However, generally speaking, arguments concerning 
the dissipation rate of W are still far from sufficient. 

Also in the context of the turbulence dynamo, the transport equation for the 
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turbulent cross helicity has been considered, where an algebraic model for the cross- 
helicity dissipation rate has been proposed.^ [3] Sur and Brandenburg wrote down 
an evolution equation for the cross-helicity effect, and argued the cross-helicity 
destruction in the context of quenching mechanism. (22J 

In order to examine the evolution of the turbulent cross helicity W, it is indis- 
pensable to properly estimate the dissipation rate of W, Ew, as well as the cross- 
helicity production rate P\y- We address this problem; modeling the cross-helicity 
dissipation rate on the basis of a statistical analytical theory. 

Spacecraft observations of solar-wind turbulence have revealed detailed informa- 
tion on the large-scale behavior of turbulent statistical quantities, which includes 
the radial evolution of the cross helicity both in the low- and high-speed wind re- 
gions. Comparison of the satellite observations with the numerical simulation with 
the aid of a turbulence model provides a good test for the cross-helicity dissipation 
models. In this work, we will delve into the problem of cross-helicity dissipation 
modeling by using such comparisons. 

The organization of this paper is as follows. After briefly showing the evolution 
equation of the cross helicity in Section [2j we present the exact equation of the 
cross-helicity dissipation rate Ew in Section |3j Then, we consider two candidates 
for the model of Ew in Section [4j one is the algebraic model and another is a 
transport-equation model. These expressions are systematically derived with the 
aid of a statistical analytical theory of inhomogeneous turbulence. Some features 
of the adopted turbulence model, which is an expansion of the hydrodynamic k — e- 
type one-point turbulence model in the engineering field, are noted in Section [5] 
An application of the model to the solar wind is presented in Section [6} A brief 
summary is given in Section [7j 



2. Equation for the cross helicity 

The fluctuation velocity and magnetic field, u' and b', in incompressible magneto- 
hydrodynamic (MHD) flow are governed by 

d + U • V ) u' = - (u • V) U + (B • V) b' + (W ■ V) B 



at 

- (u' • V) u' - (b' • V) b' + v • n - Vp' M + i/VV + f,(l) 

(^ + U • v) b' = - (u' • V) B + (B • V) u' + (b' • V) U 

- (u' • V) b' - (b' • V) u' - V x E M + AV 2 b' (2) 

with the solenoidal conditions 

V • u' = V • b' = (3) 

[y: kinematic viscosity, A: magnetic diffusivity) . The primed quantities denote the 
deviations from the ensemble average (• • • ) as 

(p = 7p + cp', Tp= (tf) (4) 
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with 

ip = (p,u,u,h,i,p,p M ,{), (5a) 

Jp = (p,XJ,n,B,J,P,P M ,F), (5b) 

V 1 = {ft ,d ,u' ,b' ,$' ,p' ,tf u ,f') . (5c) 

Here, p is the density, u(= V x u) the vorticity, j(= V x b) the electric-current 
density, p the gas pressure, Pm{= P + b 2 /2) the MHD pressure, and f the external 
force. Note that the magnetic field etc. are measured in the Alfven-speed unit. 
They are related to the ones measured in the original unit (asterisked) as 

i "* • J* e* p* , . 

b -(^) I ^ , ''IpJ^W 2 ' e "W) I ^ , p "7 w 

(/Uo: magnetic permeability). 

The Reynolds stress TZ and the turbulent electromotive force Em represent the 
turbulence effects on the mean field. They are defined by 

n * = (u ,a u'P - b' a b' p ) , (7) 

E M = <u x b') , (8) 

respectively. 

In order to describe the properties of turbulence, we consider several turbulent 
quadratic statistical quantities. Among them, the turbulent MHD energy K and 
the turbulent cross helicity W are defined by 

K^(u' 2 + b' 2 >, (9) 



W = (u'-h'). (10) 

From Equations M and pL the evolution equations of K and W are rightly 
obtained as 

tM! +u - v ) g = pg -^ +t ° (n) 

with G = (K,W). Here, Pq, £q, and Tq are the production, dissipation, and 
transport rates of the turbulent statistical quantity G. They are defined by 

r)JJ a 
P K = -K ab —-n M -J, (12a) 

ox 



du' a \ 2 \ . I fdV' 



« = "<!&? +A U? 1 >" (12b) 
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T K = B- VW - V 



+ p' M )u'+(u'.b')b^ + (f'.u'> 



(12c) 



dB a 
Pw = — H a -—£■ — Em • fi, 



(13a) 



e w = (v + A) 



du' a db' a 
dx b dx b I ' 



(13b) 



■w 



B-VK-V ■ 



u 



b')u' 



u 



n 



+ b 



n 



p'm 



W) + (f'-b') 



(13c) 



As was already mentioned in Section [T] and Appendix |Aj the production rates 
Pk and Pw represent the supply of K and W, respectively, due to turbulence 
cascade. We see from the first term of Equation ( |12a ) that the turbulent energy 
is sustained by the mean velocity shear dll a /dx b . The second or J-related term 



of Equation (12a) is related to the Joule dissipation of MHD turbulence. Equa- 



tion (13a) shows that inhomogeneities of the mean fields coupled with the fluctua- 



tion correlations are essential for Pw . One is the shear of the large-scale magnetic 
field, dB a /dx b , coupled with the Reynolds stress JZ: 



P\VK 



-n 



ab 



dB a 



0:1 



(14) 



The other is the large-scale vortical motion fl coupled with the turbulent electro- 
motive force Em: 



Pwe 



Em • O. 



(15) 



Finite positive (or negative) values of P\vr and Pwe infer the generation of positive 
(or negative) turbulence cross helicity W by way of cross-helicity cascade. The role 
of the cross-helicity production have been argued mainly in the context of turbulent 
dynamo and turbulent transport suppression. 

Other important mechanisms that possibly supply the cross helicity to turbulence 



come from the first and last terms of Equation (13c): 



T WK = B • VK, 



(16) 



T WF = (f • b') 



(17) 



Equation (16) shows that the inhomogeneity of K along the large-scale magnetic 



field B may contribute to the supply of the cross helicity. We often meet such 
situations in astrophysical phenomena where inhomogeneous turbulent plasma is 
threaded through by the ambient magnetic fields. 
On the other hand, Equation (17) represents the turbulent cross-helicity gen- 



eration due to the coupling of the external forcing and magnetic fluctuation. As 
we see in the evolution equation of the mean cross helicity U • B [Equation ( A9 ) 
in Appendix Kj , the coupling of the (mean) external forcing with the large-scale 



August 1, 2011 0:36 Journal of Turbulence n'yokoi'JOT-2010-0067'proof 

6 N. Yokoi 



magnetic field provides the mean cross helicity. Actually, such cross-helicity gen- 
eration due to an external forcing plays a crucial role in the dynamo action or 
magnetic-field generation mechanism in a generalized Arnold-Beltrami-Childress 
flow called the Archonitis flow. [22] In this sense, analysis of the forcing effect on 
the evolution of the cross helicity is very important. It may affect the expression 
of the cross-helicity flux, just as the body- force effects such as the buoyancy, frame 
rotation, etc. change the scalar transport. However, this topic is beyond the scope 
of this paper. 

Unlike the production rates Pk [Equation ( |12a )] and P\y [Equation (13a)] and 
the first or B-related terms of the transport rates [Equations (12c) and (13c)], the 



definitions of the dissipation rates, e [Equation (12b)] and ew [Equation (13b)], 



do not directly contain any mean fields, but are expressed by the combination of 
the molecular viscosity v and magnetic diffusivity A with the inhomogeneity of the 
fluctuation fields. So, we need some specific treatments on them. We consider the 
cross-helicity dissipation rate in the following sections. 



3. Equation for the cross-helicity dissipation rate 

From the equations of the fluctuation velocity and magnetic field, we construct the 



exact equation for the dissipation rate of the cross helicity, ew [Equation (13b)], 



MS 



Dew __ 
Dt ~ 

= (u + X) 
+(^ + A) 
+(^ + A) 
+(^ + A) 
-(iz + A) 
-(iz + A) 
-(iz + A) 
+(^ + A) 



0_ 

Ft 

du' a db' b 



+ U • V ] e w 

Qh ta Q u tb 



dU a 



dx c dx c dx c dx c / dx b 



rlii' a 

OU lib 

— — u 

ox c 

du' b d 2 u' b 



, b db" 



d 2 U a 



dx c / dx b dx c 



+ 



db' b d 2 b' 



HUb 



dx a dx a dx c dx a dx a dx c 
du' c du' c db' c db' c \ dB a 



dx a dx b 
du' a du' b 



Q x a Q x b J Q x b 

db' a db' b \ dB a 



dx c dx c dx c dx c / dx b 



du' a 

0:1 



-U 



ib 



db' c 
dx c 



lib 



d 2 B° 

dx b dz 



B ( 



du lb du lc db' b \ 
~dx^~dx^~dx^ J ~ ^ V ' 



du' b db' c db' b 
dx a dx a dx c 





is + \)( (u' c ± b' c ) 



+(^ + A) 



_d_ 
dx c 



dx c 

d 

dx a 



+ (^ + A) 

du' b db' b 

dx a dx a 



u' b ± b' b 



db' b du' c du' b 
dx a dx a dx c 

db lb db' c db' 1 
dx a dx a dx c 
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i-MS^W + A)^,,, 



dx a dx a dx b 



dx c dx c 



(i/ + A) 
(i/ + A) 



cte c 

2 



du' b d 2 b' b 
dx a dx a dx c 

Q2 u ,b Q2 b ,b 



+ A 



db' b d 2 -' b 



u 



dx a dx a dx c 



dx a dx c dx a dx c 



(18) 



(both upper or both lower signs should be chosen in the double signs). This equa- 
tion, lacking the connection with a conservative law, has a considerably complicated 
structure. This is in sharp contrast to the K and W equations [Equation ([XT])]. 

In the hydrodynamic case with an electrically non-conducting fluid, the energy 
dissipation rate 



la 



du' a du 
dx b dx b 



(19) 



as well as the turbulent energy K u = (u /2 ) /2, plays a central role in turbulence 
modeling. [23] From the equation of the velocity fluctuation, we write the e equation 
exactly as 



De 
Dt 



d_ 

Of 



+ U-V e 



-2v 



lb 



du' b du' c du' 
dx a dx a dx c 



} 'du' a du' b du' a du' b 
dx c dx c dx c dx c 

du' a \ d 2 U a 




-2v ( W 



,'b 



dx c 




I dp 1 du' a 
\dx b dx b 



+v 



8 2 t 



dx c dx c 



(20) 



The mathematical structure of this equation is also complicated because of the lack 
of the connection with a conservative law. 

The energy dissipation e itself is dominant at small scales. Using this, the length 
and velocity scales are estimated as 



Ixl ~ ^/V 1 / 4 , 



I /i 1/4 1/4 

u ~ v ' e ' , 



(21) 



respectively. Using Equation (21), we can estimate each term in Equation ( |20[ ). In 
the flow at high Reynolds number (y — > 0), two terms behaving as 0(y~ l l l e'l 2 ) 
are dominant, and these two terms should balance each other [3 [2"^] 



-2v 



ib 



du' a du' a du' 
dx b dx c dx c 



Zn.'a 



d 2 w 



dx b da 



(22) 
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In other word, in the modeling of the e equation, it is of crucial importance to 



properly estimate the two terms in Equation ( 22 ) 



In the hydrodynamic turbulence modeling, an empirical model equation for e: 

^ s (! +U . V ) e = Cll lft- C , 2 i e + v(^V £ ) (23) 

was proposed and has been widely accepted as useful. 

Here, C t i, C e 2, and a e are model constants. The values of these constants have 
been optimized through various applications of the k — e model. Usually, the values 
of 

C el = 1.4, C e2 = 1.9, a e = 1.0 (24) 

are adopted. [28J 

Equation ( |23| ) is based on an empirical inference that the dissipation should 
be larger where the turbulence is larger, and is derived by dimensional analysis. 



However, it is also pointed out that the system of constants given by Equation (24) 
is not a unique combination for the model constants because e equation can take 
an infinite number of self-similar states. [25] 



A similar argument using Equation (21 ) can be applied to the exact equation for 



the dissipation rate of the turbulent cross helicity [Equation (18)]. The difference 
between the molecular viscosity v and the magnetic diffusivity A is expressed by 
the magnetic Prandtl number: 

Pm = v/X. (25) 



If we estimate each term in Equation ( 18 ) under the simplest possible condition of 



Pm = 1, we have a dominant balance between the terms expressed by 

x . /du' b du' c db' b \ 1/2 , /2 , , 



(" + >■? ( f"' b fK ) ~ v-V's** (26b) 

v ' \ £i~~a i2i~,r q_a o_/~ / \ / 



% dx a dx c dx a dx 
in the Ew equation at the high Reynolds number flow. 

4. Models for the dissipation of the turbulent cross helicity 

4.1. Algebraic model 

As was mentioned in the previous section, the equation governing the dissipation 
rate of W, E\y, is very complicated. In such a situation, the simplest possible model 
for e\y is the algebraic approximation as follows. 

Using the turbulent MHD energy K and its dissipation rate e, we construct a 
characteristic time scale of turbulence as 

r = K/e. (27) 
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Other choices of time scale are possible. A large-scale magnetic field may alter 
the characteristics of turbulence. In the case of MHD turbulence, the Alfven time 
associated with the magnetic field may modulate time scale of turbulence. This 
point will be referred to later at the end of Section [5] 

With the aid of the time scale of Equation ( |27[ ), the dissipation rate of W can 
be modeled as 

W e 

e w = C w — = C w -W, (28) 

r K 

where Cw is the model constant. Namely, we consider the dissipation of the tur- 
bulent cross helicity is proportional to the turbulent cross helicity divided by the 
time scale. 

It is worth noting a mathematical constraint on the cross helicity. Namely, the 
magnitude of the turbulent cross helicity W is bounded by the magnitude of the 
turbulent MHD energy K as 

K ~ <u'2 + b'2>/2- l yj 



This relation constrains the value of e\y- Actually, from Equation (11), the turbu- 



lent cross helicity scaled by the turbulent MHD energy, W/K, is subject to 
D W W ( 1 DW 1 DK" 



DtK K \W Dt K Dt 

W ( 1 „ 1 \ W ( ] 1 

-Pw ~ T?Pk - -77 [ 7Tr e W ~ 77 £ 



K \W K ) K \W K 

+ 1 (w T * ~ » ■ < 30 > 

Equation ( |30[ ) is a very general expression for the evolution equation of the scaled 
cross helicity, W/K, which should be satisfied with in any situations of turbulence 
flow. If we consider homogeneous turbulence, where the spatial variations of mean 
quantities vanish, we have neither production nor transport rate (Pq = Tq = 0). 



In such a case, Equation (30) is reduced to 



d W ( 1 1 \W 

dtK=-\w eW -K £ )K- (31) 

This may give a constraint on the values of the dissipation rates of the turbulent 



MHD energy and cross helicity, e and Ew- At least, Equation (31 ) provides us with 



a constraint for the turbulence modeling of e and ew as we see in the following. 



If the coefficient of W/K or parenthesized quantity in Equation (31) does not 
depend on W/K, which occurs naturally when we adopt simple algebraic models 
for e and £w as below, in order for the magnitude of the scaled cross helicity to 
be bounded in time evolution, the parenthesized quantity should be positive. Then 
we have inequality 

|*| >m. (32, 



£ K 

Provided that the difference between the time scales of W and K is independent of 
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W/K, this is a necessary condition for £\y in relation to e arising from the upper- 



boundedness of W [Equation (29)]. Actually, if we adopt simple algebraic models 
for e and £\y such as 

K W 

e = — , £w = Cw — i (33) 

r r 



the parenthesized quantity in Equation (31) does not depend on W/K. Then we 



have inequality (32), which gives a constraint on the model constant C\y as will 



be seen later in Equation (88) in Section 6.2. 

In the context of homogeneous MHD turbulence, since the pioneering work by 
Dobrowolny, Mageney &; Veltri[T0l lllj . it had been considered that the initially 
prevailing cross helicity grows with time due to the nonlinear interaction towards 
the +1 (or —1) value of the scaled cross helicity (Also see |12|I13|). However, fur- 
ther investigation by Ting et al.[26] and Stribling & Matthaeus[27] showed that 
MHD turbulence has diverse possibilities for long time evolution, including either 
growth or reduction of W/K. In contrast to such homogeneous MHD turbulence, 
in inhomogeneous MHD turbulence with non- vanishing mean velocity shear, we 
have production, dissipation, and transport mechanisms of the turbulent cross he- 
licity arising from or related to the inhomogeneity of mean fields. In the present 
work, we explore these mechanisms intrinsic to the inhomogeneities of MHD tur- 
bulence. There have been several works where the effects of mean velocity shear 
are incorporated into the evolution equations of turbulent quantities including the 
turbulent cross helicity. In this sense, this work is in the same line with |18H20| . 
and also |28H30] . However, in the present work, special emphasis is placed on the 
theoretical derivation of the turbulent cross-helicity dissipation-rate equation as 
will be shown in the following sections. 

Here, we had better remark on the notation of the scaled cross helicities. The tur- 



bulent cross helicity scaled by the turbulent MHD energy, W/K [Equation (29)], is 
"dynamically" important in the context of turbulent dynamo etc. There is another 
scaled cross helicity 

- <»'"'> (34) 

V(u' 2 )(b' 2 ) 

that is "kinematically" or geometrically important, since it represents the degree 
of alignment of the velocity and magnetic-field vectors. In the solar-wind turbu- 
lence community, the former, W/K, is often referred to as u ac" ■ However, in the 
"general" turbulence community, this is not the case; Some use "pc" f° r this quan- 
tity, others do "p", and so on. Hence, in this paper, we confine ourselves to just 
mentioning these notation conventions. 

Interestingly, the directional alignment expressed by T has shown to be robust 
rather than the behaviour of W/.K".[3~I] This tendency is related to the fact that, 
unlike T, W/K may dynamically change its value depending on how much differ- 
ently the turbulent MHD energy and cross helicity, K and W, are influenced by 



the large-scale shears. This point will be referred to later in Sec. 6.4 



4.2. Model equation for the cross-helicity dissipation rate (ew equation) 

The equation of the K dissipation rate, e equation, can be written as 

|? + (U • V) £ = C El ^P K -C £2 ^£ + V- (^ Ve) . (35) 
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This is a direct expansion of the hydro dynamic (HD) or non-MHD turbulence 



energy dissipation e equation (23). The model constants in Equation (35) are same 



as the ones appearing in Equation ( 23 ) . This is a consequence of the requirement 
that the e equation for MHD turbulence should be reduced to the e equation for 
the HD in the limit of the vanishing magnetic field (B = 0, b' = 0). 



We see from Equation (11) that the equations of W and K are written in a 
similar form. So, it is natural to consider that the equation of the dissipation rate 
of W, Ew, can be expressed in a form similar to the equation of the K dissipation 
rate e. Then we may consider the equation for ew as 

—f^~ = C W1 —P W - C W2l? e w + V • ( — —Ve w ) , (36) 

where Cyyi, Cw2, and a £ w are model constants. 

In contrast to the MHD energy, the cross helicity, defined by the correlation 
between the velocity and magnetic field, vanishes in the limit of the vanishing 
magnetic field. As this result, as far as the ew equation is concerned, we can not 
make use of the knowledge accumulated in the history of HD turbulence mod- 
eling. To say nothing of the model constants Cwi, Cw2, and a £ w appearing in 



Equation ( 36 ) , we have to consider the structure of ew equation itself from more 
fundamental basis. For this purpose, in the following we will construct the equation 
for the cross- helicity dissipation rate on the analytical theoretical basis. 

From the theoretical analysis of inhomogeneous MHD turbulence, the cross- 
helicity density defined by W = (u' • b') is expressed as 

W = 2I {Q ub } -iJgs,^ {Q ub + Q bu )\ , (37) 

where we have used abbreviated forms of spectral and time integrals: 

I n {A} = J A(k, x; r, r, t)k 2n dk, (38a) 

I n {A,B}= dkk 2n dTiA(k,yL;T,T h t)B(k,xL;T,Ti,t). (38b) 

J J -co 



In Equation (37), G§ is the Green's function, and Q ub and Q bu are the correlation 



functions related to the basic or lowest-order fields u' B and b' B : 



;u B -b B ) = 2 J Q ub (k;r,r)dk. (39) 



For the derivation of Equation (37) and higher-order expression for W, see Ap- 
pendix |B} Suggestions from the higher-order expressions for the cross- helicity dis- 
sipation model are also presented in Appendix [Cj 

We assume that the correlation function and the Green's function in the inertial 
range are expressed as 

Q ub (k,x;T,T',t) = a w (k,x;t)exp [-w w (k,x; t)\r - t'\] , (40) 

Gs(k,x;T,T',t) = e(T-T')e W [-Lu s (k, X] t)(T-T')], (41) 
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where 9{t) is the Heaviside's step function that is 1 and for r > and r < 0, 
respectively. Here, aw is the power spectra of the turbulent cross helicity, and uiyy 
and wg represent the frequencies or time scales of fluctuations. As for the spectra 
in the inertial range, we assume 



a w (k,x;t) = awo£ 1 ^ 3 £w('X-;t)k n ' 3 , 
and for the time scales, 

ws(fc,x;t) = Lu S0 e 1/3 k 2 / 3 = t s \ 



(42) 



(43) 



ui w (k, x; t) = u W Q£ w Z k 2/3 



•\Y 



(44) 



{(Two, wso, and ojwo are numerical factors). Equation (42) arises from the assump- 
tion that the spectrum of the cross helicity is determined by the scale, energy and 
cross-helicity transfer rates; k, e and £w- 



Using Equations (40) and (41), W is calculated as 



W = 2 f dk ^ 1 / 3 (x;t)e V i/(x;t)^ 11/3 

1 Daw cw Dujw 



2 / dk 



UJS+Uw Dt (uj s + ujw) 2 Dt 



(45) 



From the inertial-range forms (42)- (44), this can be rewritten as 



W = 4 • 2ire- 1/3 e w I dk k~ 5/3 
J\k\>k c 



-4-27T 



&wo 



W S w£sw ^|k|>fc 



dk k A ' 3 



D 
Dt 



e 1/3 {pc;t)e w {x;t)k 



-11/3 



. , n VWO^WO -1/3 f r, , -3 D 

+4-2tt — — e l ' e w / dk k d — - 



-1/3^. 

--w 



K-t)k 2 ' 3 



(46) 



{kQ\ cut-off wave number). Here, we have introduced a synthesized time scale rsw 
defined by 



1 1 1 

+ 



TSW T S TW 



u S0 e 1/3 + uj W0 s% 3 ) k 2 ' 3 = cj sw eli 3 k 2 / 3 . (47) 



1/3 

We should note that w sw and e sw will appear only in the combination of o; sw e s 4 ; 
each of uj sw and e sw has no definite meaning. For simplicity of notation, hereafter 
we denotes 



1/3 



At \ - ^wo^w 
Aw{US0,UW0) = Y/3 



1/3 
UWO^w 



T~SW 



u sw £ s i a; S oe 1/3 + u W o£ w Tw 



(48) 



In Equation (46), the lower bound of the spectral integral region, kc, is directly 
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connected to the largest eddy size of turbulent motions as 

£ c = 27r/k c , 



13 



(49) 



which dominantly contributes to and determines the turbulent MHD energy and 
cross helicity. Strictly speaking, the correlation length of the cross helicity, lw, can 
be different from that of the MHD energy, Ik- However, if we consider the fact that 
the characteristic lengths of turbulence represent the scales of the largest turbulent 
motion corresponding to the scales with largest magnitudes of the turbulent MHD- 
energy and cross-helicity spectral densities, we see that in order for the length 
scales of the MHD energy and cross helicity to be considerably different from each 
other, the spectral distributions of the MHD energy and cross helicity should be 
entirely different from each other. This is not the case, for instance, in the case 
of the solar-wind turbulence. Actually, in most cases of interests, these two length 
scale are similar to each other. In this sense, we can regard the characteristic 
lengths of the energy and cross helicity are approximately the same, and we denote 
£ w ~ £k(= £c) hereafter. In this respect, extensive discussions of correlation length 
scales as well as time scales found in Matthaeus et al.[32j and Hossain et al.[33J 
are very important. 
With this point in mind, we perform the Lagrange derivatives and calculate the 



spectral integrals in Equation (46). Denoting the scaled wave number 



s = k/kc, 



(50) 



we have 



W 



4 • (2^) 1 / 3 cr H /oe" 1/3 ewA;c /3 / ds s~^ 3 



1 2tt aW0 k 7/3 / dss- 7 ' 3 - 
4 2vr k c dss 

W sw £sw 



«>1 

D 



S>1 



.-1/3 



e i/d {x;t)e w {x;t)k c 



11/3 



. „ &W0U W q _ 1/3 2 / 

+4 ' 27T : I73^ e £ wk c J dss 



-7/3 



s>l 



D 
Dt 



-1/3 



£ w ( x ; ^)^c 



2/3' 



Using Equation (|49|), we calculate Equation (|5 1 p to obtain 
W 



6-(27r) 1 / 3 cT W oe _1/3 £VK4 /3 



(51) 



1 



&W0 



(27r)V3 



1/3 

Wsw^sw 



WewtH* 



IDs 

~e~Dt 



[3 — Aw(coso,^wo)] 



1 Dew 
Ew Dt 



1 D£ c ] 

+ [11 -2A w (aj so ,u}wo)] -f~Dt I ' 



(52) 



In hydrodynamic turbulence modeling, the turbulent energy K u (= (u /2 )/2), its 
dissipation rate e, and the correlation or integral length £q are equivalently impor- 
tant. For the purpose of closing the system of model equations, we can select any 
combination of K u , e, and lc- This is known as the transferability of the model 
with respect to K u , e, and £q. In order to satisfy this transferability requirement, 
these quantities should be connected with each other in an algebraic relation: 



K u = C K e 2 lH 2 '\ 



(53) 
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This transferability requirement gives a theoretical foundation of the energy- 



dissipation-rate equation (23).|34j 
Relation (53) itself has been argued for a long time. As for the theoretical and 



experimental basis for this relation, the reader is referred to classical papers [351136] 
and also to Batchelor's book [37] and works cited therein. As for the simple physical 
arguments for this relation, see [3"8"H4"U] . Here, we just point out that this relation 
is easily obtained if we assume the local equilibrium of turbulence; Pku — e with 
the notion of the mixing length. Then we can estimate the dissipation rate as 



la„,lb 



Pku = ~{u' a u 



dU b 
dx a 



~ uu- 



(54) 



(u: characteristic intensity of turbulence, £: mixing length). This is equivalent to 



Equation (53). 



In the context of theoretical derivation of e equation, it is the algebraic property 



of Equation (53) that is much more important. If the relation is not algebraic, we 



have no transferability among any combination of K u , e, and £q at all. 
We expand the transferability requirement to the model equation related to the 



cross helicity. For this purpose, we solve Equation (52) concerning £q in a pertur- 
bational manner. At the lowest-order, we have 



W = 6(2Tr) 1 / 3 a W0 e- 1/3 e w l 2 ( / 3 



(55) 



or 



£ C = 6- 3 / 2 (27T)-VV w f ^^3/2^3/2. (56) 

The transferability of model requires an algebraic relation among W, e, ew, and 



£c- Equations (55) and (56) correspond to such a relation. 



Using Equation (|56|), we change expression (|52|) based on e, ew, and £q into the 

have 

(57) 



one based on e, ew-, and W. We require Equation (55) or (56). As a result, we have 
Dew 



, e w De e w DW 

— — - Ci{uso,uj wo ) — — — + o 2 (wso,wwoJtt^-=— 



with 



n i % 13 — -Avp(wso,wwo) n i \ 33 - 6A w (ojso,^wo) 

<^l{USO,Uwo) = ™ ^Tl 7~. ~ \i ^2(^30,^1^0; ~" 



39 - 8A w (ujso,^>wo) 



39 — 8Aw(u}so,uwo) 



(58) 

If the cross-helicity time scale can be treated as similar to the time scale associated 
with the Green's function: 



1/3 1/3 

wsoe ' - u W q£ w , 



Awi^so, wwo) [Equation (48)] can be estimated as 

Awi^so, uwo) — 1/2. 
Then we have model coefficients 

Ci(u> So ,u>wo) - ~ - 0.34, C 2 (uj S o,ujwo) - ■= - 0.86. 
35 7 



(59) 



(60) 



(61) 
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If we substitute equations of e and W into Equation (57), we finally obtain 



Dew r, £\V D n £W . r, £\V D ^ £ W , coA 

with 

Ce w l = ClCel) Ce w 2 = C\C e \, C £w 3 = C £w 4 = C2 (63) 



[C e i and C e 2 are given as Equation (24)] 



From Equation (55), we shortly have 



^ = 6- 1 (2vr)- 1 /V^ 1 ^ (64) 

with 

T = e- 1/3 e 2 J 3 . (65) 



This corresponds to the simplest possible algebraic model [Equation (28)] for the 



cross-helicity dissipation rate. In other word, the simplest possible algebraic model 



for the cross-helicity dissipation corresponds to the scaling of Equation (42). 

In the presence of a mean magnetic field, turbulence deviates from the isotropic 
state and the resulting anisotropic effects should be taken into account (See e.g., 
|41j). As is shown in Appendix ICj in the present formalism, we can take the effects 
of mean magnetic field and large-scale rotation into consideration, and obtain the 
the cross-helicity dissipation equation with such effects incorporated. Related to 
this point, we should make the following remarks. In this work, we assume an 



isotropic form Equation (42) for the basic or lowest-order fields, u' B and h' B , but 
the deviations from homogeneous isotropic state are taken into consideration from 
the higher-order contributions related to the inhomogeneity of mean velocity and 
the presence of the mean magnetic and vorticity fields. In this sense, the turbulence 
fields themselves are neither isotropic nor homogeneous in this work. However, 
this point never denies the importance of incorporating the anisotropy effects into 
the basic fields. This is a very interesting topic in the theoretical formulation of 
inhomogeneous turbulence, and will be reported in the forthcoming paper. 



5. A model for MHD turbulence 

One of the prominent features of turbulence is its wide ranges of scales. We have 
continuous scales of motions ranging from the energy-containing scale, in which 
energy is injected to turbulence through the inhomogeneity of mean fields, to the 
dissipation or Kolmogorov scale where the dissipation plays a dominant role. Owing 
to this breadth of scales, it is impossible in the foreseeable future to simultaneously 
solve all the scales of turbulence at the high Reynolds number encountered in the 
real-world flow of interests. In such a situation, the notion of turbulence modeling 
provides a useful tool for analyzing the real-world turbulence phenomena. In the 
turbulence modeling, small-scale motions are modeled, and turbulence effects are 
incorporated into the analysis of the large-scale or mean motions. 
In an MHD fluid, the mean fields obey 
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Equation of continuity: 



dp 
di 



+ V • (pU) = 0, 



(66) 



Momentum equation: 



dt 



+ V • (UU - BB) = - VP M -V-n + vV 2 XJ, 



(67) 



Magnetic induction equation: 

V x (U x B) + V x E M + AV 2 B, 



<9B 

~dt 



(68) 



Solenoidal condition for the magnetic field: 

1 



V B 



(B-V)lnp. 



(69) 



Here, p is the mean density. For the sake of simplicity, the fluctuation part of the 
density has been neglected (// = 0). Measured in the Alfven-speed unit [Equa- 
tion Q] , the solenoidal condition of the original mean magnetic field: 



V • B* = 



(70) 



is expressed as Equation ( 69 ) . 



Equations (67) and (68) show that the Reynolds stress 1Z. and the turbulent 
electromotive force Em defined by Equations ffh and ^ , represent the turbulence 
effects in the mean-field equations. Estimating these correlations is of central im- 
portance in the study of inhomogeneous turbulence. In order to close the system 
of equations, they should be modeled in terms of the mean-field quantities. Here, 
we adopt 



n 



a/3 



K K 5 ap - u K S ap + v M M 



ap 



a a/3 



(71) 



E M = aB - /3V x B + 7 V x U, 



(72) 



where S and A4 are the strain rates of the mean velocity and the mean magnetic 
field: 



**-w+£-\r*-* 



(73a) 



*-- %^- 2 r— 



(73b) 



We should note that the transport coefficients appearing in Equations ( 71 ) and 



(72), z/jc, z^m, o>, f3, and 7, are not adjustable parameters, but should be determined 



self-consistently through the dynamic properties of turbulence. They in general 
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depend on location and time. In this work, we adopt model expressions for them 

as 



5 r ic 5 

-UK = CprK, 7 = -i 



a = C a rH, p = -v K = C p tK, 7 = -v M = C^tW, (74) 



where 

H=(-u'-u;' + W-i'), (75) 

and r denotes the time scale of turbulence. The model constants C a , Cp, and C 1 
are estimated as [5] 

C a ~ 0.02, Cfi ~ 0.05, C 7 ~ 0.04. (76) 

In the Reynolds- averaged turbulence model, properties of turbulence are repre- 
sented by some statistical quantities. We adopt the turbulent MHD energy K, its 
dissipation rate e, the turbulent cross helicity W, the cross-helicity dissipation rate, 
ew, and the turbulent residual energy Kr as the turbulence statistical quantities. 



The definition of them are given by K [Equation (|9j)], W [Equation (10)], 

^R - \ <u' 2 - b' 2 > , (77) 

and the dissipation rates of the turbulent MHD energy and cross helicity, e [Equa- 



tion (12b)] and Ew [Equation (13b)], respectively. 



In this work, we assume the a- or turbulent residual-helicity-related term in 



Equation ( 72 ) is much smaller than the 7- or cross-helicity-related term, and neglect 
the former. So, we do not include the evolution equation of H. As for the modeling 
of H equation and inclusion of it into the turbulence model, see [5l [301 H2] • Related 
to this point, the equations of kinetic helicity in HD turbulence and in MHD 
turbulence are an interesting subject. For these topics, the reader is referred to |l3] 
for the HD case and [H] for the MHD case. 

If we take the effects of mean-density variation into account, the evolution equa- 
tions for the statistical quantities are written as 

^ = - (U • V) K - \ (3K + K R ) V • U 
at 6 

+ ^{v K S 2 - v M S : M) + /3J 2 - 7ft • J 

-e + V • (WB) + V • {v K VK) , (78) 



dW 1 

— = -(U. V )W--WV-XJ 

+- (u K S : M - is M M 2 ) +(3Q-J - 7 n 

-e w + V-(KB) + V ■ (— VW) , (79) 

\crw J 



2 
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^ = -(\J.V)K R -±(K + 3K R )V-lJ 



+ l - UK (S 2 -M 2 )-e K + V 



UK 



VK R 



(80) 



where S 2 = S ab S ab , S : A4 = S ab M ab , etc. Here z/r denotes the residual turbulent 
viscosity defined by 



VR = V K 



Kr 
K 



(81) 



We should note that Equations (78) and (79) are very similar to equations found 
in Matthaeus et al.[32j 
The dissipation rate of the turbulent energy, e, in Equation (1781 is given by Equa- 



tion (35). The dissipation rate of the turbulent cross helicity, E\y [Equation (13b)], 
is estimated by Equation ( 28 ) or alternatively, by Equation (l62|) . 



The dissipation rate of K R , £r, in Equation (80) can be expressed as [28 



^R 



c 4)k k >' 



(82) 



Here, C r and cjr are positive model constants. The large-scale behavior of K R 
depends on the choice of these constants. At this stage of modeling, we do not 
insist on fine tuning of the model constants, but roughly put them 



C r = 0.01, CTR = 1.0. 



(83) 



Equation (82) as a whole will destroy the turbulent residual energy, and re- 



turn MHD turbulence to equipartition. The first part in the parentheses of Equa- 



tion (82) represents the K R destruction due to the eddy distortion. The second 



part represents the destruction of K R due to the mean magnetic field. The latter 
corresponds to the Alfven effect, in which the presence of the mean magnetic field 
leads MHD turbulence to an equipartition state between the kinetic and magnetic 
energies. |28| H5rl47j In connection with the anisotropy and residual energy, sev- 
eral numerical studies have recently been performed in the presence of a uniform 
magnetic field. For this interesting topic, see |48H50j and works cited therein. 



Equation (82) may be reinterpreted as the modulation of MHD-turbulence time 



scale due to the mean magnetic field: 



K 

£ 



i + x 






-i 



K 

e 



(84) 



where x(= Cr) represents the synthesization ratio of time scales. Equation (83) 
infers that we should put the ratio as 



X = O(10~ 



(85) 



As for an example considering the synthesized time scale of MHD turbulence, you 
are referred to |30] and references cited therein. 
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6. An application to the solar-wind turbulence 

6.1. Solar-wind turbulence 

Solar wind is a continuous plasma flow blown away from the coronal bases to the 
solar-system space. Its origin is considered to be the violent magnetic activities on 
the solar surface. The streams from the low- or mid-latitude region whose typical 
speed is 400 km s~ are called the slow wind. On the other hand, the streams from 
the high-latitude region such as coronal hole whose typical speed is 800 km s~ are 
called the fast wind. It is known that the slow wind has a large velocity shear, while 
the fast wind has substantially no velocity shears. The influence of an explosive 
magnetic activity on the solar surface is conveyed to the magnetosphere of the Earth 
within several days, and terrestrial environments may be affected much. A direct 
exposure to a high-energy particle convected by the solar wind is quite harmful to 
astronauts and electric devices on the satellite. So it is highly desirable to predict 
the behavior of the solar wind: A problem of the space weather (Figure nj). 




Figure 1. Solar wind and terrestrial environments. The solar wind, a high speed plasma flow of sev- 
eral hundreds km s , interacts with the terrestrial magnetosphere to induce several phenomena such as 
substorms, auroras. 



Satellite observations have revealed the statistical properties of the solar-wind 
turbulence including the velocity, magnetic field, density, etc. It has been investi- 
gated how MHD turbulence evolves in the large-scale velocity and magnetic-field 
structures. |14IfT6] According to the observations, the solar- wind turbulence shows 
a strong Alfvenicity near the Sun. Namely, there is a strong correlation between 
the velocity and magnetic-field fluctuations, and equipartition between the kinetic 
and magnetic turbulent energies is realized. At the same time, it is pointed out 
that this Alfvenicity decays as the heliocentric distance increases. |16j One of the 
main unsolved problems in this field is the spatial evolutions of the turbulent cross 
helicity W and the Alfven ratio rA-[5i| 

The magnitude of the scaled turbulent cross helicity, \W\/K, whose value is 
almost unity near the Sun, decreases as the heliocentric distance increases. In the 
region with larger mean- velocity shear, the decay rate of | W|/AT is larger. The value 
of \W\/K far from the Sun is small (0 — 0.2). On the other hand, in the region 
with smaller mean- velocity shear, the decay of \W\/K is suppressed. The value of 
\W\/K remains to be large (0.4 — 0.7) even far from the Sun. [52] The previous 
models could not reproduce this large value of the scaled W in the region with 
small or almost no mean-velocity shear. 

On the other hand, the Alfven ratio ta is defined by the ratio of the turbulent 
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kinetic and magnetic energies: 






»*A = ^«f. (86) 



This ratio is almost unity near the Sun, then decreases with the heliocentric dis- 
tance. At about 3 AU (astronomical unit: 1 AU = 1.5 x 10 11 m) from the Sun, 
r\ reaches 0.5 and thereafter remains to be the same (~ 0.5) with the heliocentric 
distance as long as the observations exist: 

r A ~ 0.5 for r > 3 AU. (87) 



Previous research could not properly explain this stationary value of 0.5. [16J 

In this work, we address these problems from the viewpoint of the turbulence 
model. In particular, we examine how the large-scale spatial evolution varies as the 
model expression for the turbulent cross- helicity dissipation rate changes. 



6.2. Numerical simulations 

We examine the evolution of the turbulent statistical quantities under the pre- 
scribed mean velocity and magnetic field. The basic set-up of the numerical simula- 
tion is the same as the previous work, so, for the details of the numerical simulation 
the reader is referred to [29U30J . For the choice of the mean fields, we fully utilize the 
current theory and modeling of the solar wind. The magnetic rotator model of the 
solar wind takes the effects of rotation and magnetic field into consideration, and 
known to be a good approximation for the mean velocity and magnetic field of the 
solar wind. |53l |53j However, its large-scale velocity and magnetic-field shears are 
weak as compared with the large-scale shears in the slow-fast-wind and magnetic- 
sector boundaries. So, the adopted mean-field profile is suitable for representing 
relatively calm flow fields within one magnetic sector. 



First, we adopt the algebraic expression [Equation (28)] for the dissipation of the 



turbulent cross helicity. In this case, Equation (32) is rewritten as 

C w > 1. 



Namely, the model constant Cw for ew in the W equation should be larger than 
the counterpart (= 1) for e in the K equation. 
In this work, we adopt 

C w = 1.4, a w = 1.0. (89) 

Here, we should note the following point. In the previous work, Cw = 1-8 was 
adopted for the algebraic model constant for £vk-|2S] As was mentioned above, the 
assumed mean fields correspond to a weak-shear case in which W is expected to 
show no considerable decay. If we consider this fact, we see that the less dissipative 



value for Cw [Equation (89)] is more suitable than the more dissipative value of 



Cw = 1.8 previously adopted. 
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6.3. Results 



6.3.1. Alfven ratio 



The simulated Alfven ratio ta is shown in Figure [2] with the comparisons with 
the observations and with the previous work. 




Figure 2. The radial evolution of the Alfven ratio r&. •: observations by Roberts et al.|15|: : the 

present model; : the present model with enhanced shear. Previous work, : Zhou and Matthaeus[17|: 

— — — : Tu and Marsch 55 . For much recent work for comparison, the reader is referred to Breech et 
al.l20l 



The evolution of Alfven ratio ta does not depend so much on the value of Cw- In 
the case of higher velocity shear (the dotted line in Figure ^1) , the large scale evolu- 
tion of ta shows a more stationary behavior in space, which is in better agreement 
with the observations. As was already mentioned, if we made the fine tuning of C r 
value, the agreement of the simulation results and observations becomes better. In 
this work, however, we do not do such fine tuning, and made only a rough estimate 



of the constants as in Equation ( 83 ) 



In Equation ( 80 ) , the turbulent residual energy equation contains a term related 
to the expansion of the large-scale flow through V • U: 



Pkri 



6 



(K + 3K R )V • U. 



(90) 



In the solar- wind turbulence, the dominance of Alfven- wave effects would give an 
equipartition of the kinetic and magnetic energies near the Sun. Therefore, we have 
K R = for the inner boundary near the Sun. Since the turbulent MHD energy is 
always positive (K > 0), we have negative contribution from Pkri near the Sun: 



Pkri 



-KV -U <0, 



(91) 



which produces a negative residual energy {K R < 0), i.e., magnetic dominance 
there. Due to this negative production, the magnitude of negative K R increases 
until its value reaches 



K + 3K R = 0. 



(92) 



If the magnitude of K R becomes larger than K/3, the sign of the production Pkri 
reverses to positive, then the production of negative K R is stopped. This infers 
that K R = —K/3 is an equilibrium state of MHD turbulence in the expanding 
solar-wind geometry with equipartition at the inner boundary (Figure pi). 
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Figure 3. Expansion effect and equilibrium of the turbulent residual energy. (Left) Flow expansion. (Right) 
Schematically depicted equilibrium state of the turbulent residual energy Kr . From the equipartition state 
(Kji = 0) near the Sun, Kj^ starts decreasing towards the equilibrium value of K^i = —K/3. 

We should note the fact that the Alfven ratio of 0.5 corresponds to the K-&/K 
of —1/3 as 
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(93) 



Equation ( 92 ) is suggestive for the explanation of the observed value of the Alfven 
ratio r\ ~ 0.5. In the expanding solar-wind flow, due to the negative K-& produc- 
tion by Prri, a state of magnetic dominance is realized. However, the magnitude 
of -Kr can not exceed K/3, since the negative production associated with a flow 
expansion vanishes when K-r, reaches —K/3. Equation (92) reflects an equilib- 



rium state of MHD turbulence in this geometry. Once this state has been realized, 
the flow expansion or compressibility makes no contributions to the Alfven-ratio 
evolution, and the turbulence evolves according to the balance of the other evolu- 
tion mechanisms intrinsic to incompressible MHD turbulence. In other words, the 
value of Alfven ratio is determined by the expansion effect with the inner bound- 
ary condition. However, once turbulence has fell into the state of r& — 0.5, the 
compressibility effects disappear, and turbulence shows a state stationary in space 
according to the incompressible balances. 

6.3.2. Cross helicity 

The simulated cross helicity W is shown in Figure [4] with the comparisons with 
the observations and with the previous work. 
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Figure 4. The radial evolution of the scaled turbulent cross helicity W . o: observations by Roberts et 

a.1 . 1 1 5 1 ; : the present model; : the present model with enhanced shear. Previous work, : Zhou 

and Matthaeus|17|: — — — : Tu and Marsch 55 . For much recent work for comparison, the reader is 
referred to Breech et aI. I20| 

The radial evolution of the cross helicity in the higher velocity shear shows bet- 
ter agreement with the observations. The mean fields adopted in this work corre- 
spond to a weak-shear case, while the observations were performed in situations 
with stronger shears. Then, the tendency that the simulation results in the high- 
shear case are in better agreement is preferable. This agrees with the observational 
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findings that in the fast solar wind, which has a small velocity shear, the scaled 
turbulent cross helicity \W\/K is kept large without decaying. In the simulation 
with a weak shear, \W\/K is about 0.5 at 3 AU and about 0.3 at 10 AU, which are 
kept relatively large values. 
From this agreement of the numerical simulation with the observations, we see 



that the algebraic- model expression for e\y [Equation (28)] is good enough, as far 
as the application to the solar wind is concerned. 



6.4. Cross helicity and velocity shear 

As was referred to in Section 6.1, the magnitude of the scaled turbulent cross 
helicity \W\/K decays much in the slow-wind region where the velocity shear is 
strong, while \W\/K decays less in the fast- wind region where the shear is weak.[52J 
This basic tendency can be elucidated from the viewpoint of turbulence model 
through a simple argument. 



The behavior of \W\/K is governed by Equation (30). Among the terms, we pick 
up the production-related terms which are directly connected to the mean- velocity 
shear as 



D W 
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"w 
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Pk 



w 
Ik 



(94) 



If we take the compressibility effect into consideration, the equation of K and 



W are expressed as Equations (78) and (79). The production rates added by the 



compressibility- or V • U-related terms in Equations ( 78 ) and ( 79 ) are given as 



P K = - \ (3K + K R ) V • U + l-Cp—S 1 

6 10 £ 



7 KW 
-C 7 — <S:.M, (95) 



Pw = ~\WV • U + ^C^S : M 



7 KW 

To 67 " 



e 



M 2 , 



(96) 



where use has been made of Equations ( 71 ) , ( 72 ) , and ( 74 ) with time scale r = K/e. 
We substitute Equations ( 95 ) and ( 96 ) into Equation ( 94 ) with the assumption that 



the mean velocity shear is much larger than the mean magnetic shear: 



<S 2 > IS: Ml M 2 . 



(97) 



Then, we have a contribution from the mean-velocity shear as 



D W 
~Wt~K 



1#R 
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(98) 



The first or V • U-related term in the parentheses of Equation (98) represents 



the expansion effect. In an expansion flow such as solar wind, V-U > 0. The 
residual energy Kr in the solar wind is nearly equal to zero near the Sun, and 
K-&/K = —1/3 far from the Sun [Equation (87)]. So, the contribution from the 



V • U-related term is negative in the solar wind case. Namely, the expanding flow 
coupled with A"r < (the magnetic-energy dominance) leads to a decrease in the 
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magnitude of W/K. Equation (98) can be rewritten as 



7 r K & W j < n, q ^ 

, — C« — o -— + •••, near the bun 

D IT j 10 p e K v ; 

DtK 1 / 1 „ T„K„,\iy ! ' ,!) ' 



18 



V • U + — C«— S 2 — + .... (far from the Sun) 
10 e y K 



Since both Cp and K/e are positive, in the presence of the mean-velocity shear, 
irrespective of the sign of W/K, it works for a decrease in the magnitude of W/K. 
If we scrutinize the large-scale velocity-shear effects on the production rate of 
the turbulent cross helicity W, we see that whether W may increase or decrease 
depends on the coupling of the mean velocity and magnetic-field shears [the second 



or S : A't-related term in Equation (96)]. At the same time, the turbulent MHD 



energy K always increases in t he p resence of the mean- velocity shear [the second or 



S -related term in Equation (95)]. As this result, the magnitude of the turbulent 
cross helicity scaled by the turbulent MHD energy, |W|/if, will decrease in the 
presence of the mean- velocity shear in the primary sense. The present model reflects 
this mechanism properly. So, this model is particularly promising for describing the 
evolution of the solar-wind turbulence. 



7. Conclusion 

The evolution of turbulent cross helicity (W = (u' • b')) was investigated from 
the viewpoint of generation and destruction mechanism of it. In particular, two 
possibilities of expressing the dissipation rate of W, £\y, were presented; (i) the al- 
gebraic model and (ii) the evolution equation for ey/. It was shown that both model 
expressions can be systematically derived from the zeroth- and first-order calcula- 
tions of the statistical analytical theory of inhomogeneous turbulence, respectively. 
Validity of the model expressions was examined with the aid of a turbulence model 
constituted by four one-point turbulent statistical quantities (the turbulent MHD 
energy K, its dissipation rate e, the turbulent cross helicity W, and the turbulent 
residual energy K-r). It was shown that, as far as the application to the solar- 
wind turbulence is concerned, the algebraic model for e\y with one model constant 
gives results plausible enough. Dependence of the cross-helicity evolution on the 
large-scale velocity structure was also discussed. In the context of the solar wind, 
it was shown that both of flow expansion and large-scale velocity shear contribute 
to decreasing the magnitude of the scaled cross helicity, which confirmed earlier 
results. [18H20] More detailed expressions for the Ey/ equation were also indicated 
from the higher-order calculation of the statistical analytical theory of inhomoge- 
neous turbulence. 
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Appendix A. Mean magnetohydrodynamic (MHD) energy and cross-helicity 
equations 



The evolution equation of the turbulent cross helicity density (u ; -b')(= W) is 



written as Equation (11) with Equation (13). Since the total amount of the cross 



helicity, J v u • hdV, as well as that of the magnetohydrodynamic (MHD) energy, 
fy(u 2 + h 2 )/2dV, is an inviscid invariant of the MHD equation, we can draw a 
clear picture for the W evolution. The production of the turbulent cross helicity 
(u ; • b') arises from a kind of cascading process in turbulence as the sink or drain 
of the mean cross helicity U • B(= VFm)- 

This is an argument similar to the one for the evolution of the turbulent kinetic 
helicity density (u' • u/) in the hydrodynamic (HD) turbulence. In the HD case, 
the total amount of the kinetic helicity f v u • udV, as well as that of the energy 
JyU 2 /2dV, is an inviscid invariant. In the HD case, such a clear picture was pre- 
sented with the evolution equations of the turbulent and mean helicity densities, 
(u'-u/) and U -agS! 

In order to see the evolution picture of the MHD energy and cross helicity clearly, 
here we write down the equations of the mean MHD energy and cross helicity 
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densities: 

K M = (U 2 + B 2 )/2, 



27 



(Al) 



W M = U-B. 



(A2) 



From mean velocity and magnetic field equations (67) and (68), the evolution 



equations of the mean MHD energy and cross-helicity densities, Km and Wm, are 
obtained as 



d_ 
dt 



+ U • V Cm = Pgm - £gm + T G 



Al 



(A3) 



with Gm = (Km,W~m)- Here, Pgm, £gm, and Tqm are the production, dissipation, 
and transport rates of the mean MHD energy and cross helicity, respectively. They 
are defined by 



KM 



+K 



ab 



dU b 
dx a 



+ Em • J = —Pr, 



(A4a) 



'dU a Y x (dB a \ 

£KM = V I -TT-T I + A I r- I 



dx b 



dx b J ' 



(A4b) 



Trm = Tkmt + Tkmb, 



(A4c) 



dB b 

P\vm = +Tl a -p—[ + Em • ^ = — P\v, 



(A5a) 



e W u = {y + X) 



dU a dB a 
dx b dx b ' 



(A5b) 



T\VM = T\VMT + T\VMB- 



(A5c) 



Here in the transport-rate terms, Tkmt and Twmt are the genuine transport rates 
of the mean MHD energy and cross helicity given by 



Tkmt = V • (-U : U + E M x B) , 



(A6) 



T W mt = V • (-B :71 + EmxU) (A7) 

[(U : H) a = U b TZ ba and (B : TZ) a = B b lZ ba }. We have dropped the v- and A-related 



terms in Equations (A6) and (A7). 



On the other hand, Tkmb in Equation (A4c) and Twmb in Equation (A5c) are 



the transport rates arising from the inhomogeneities along the mean magnetic field 
B and from external force. If we have some inhomogeneity along the mean magnetic 



August 1, 2011 0:36 Journal of Turbulence n'yokoi'JOT-2010-0067'proof 

28 N. Yokoi 

field, these terms work as a kind of generation mechanism for the large-scale MHD 
energy and cross helicity. These terms can be written as 



T™ = B • [V (U ■ B)] + U ■ (F - VP! 



u 



(A8) 



TwMB = B • 



V 



U 2 + B< 



+ F- VP 



M 



(A9) 



We should note that the production rates of the mean MHD energy and cross 
helicity, Pkm [Equation ( |A4a )] and Pwm [Equation (A5a)], are exactly the same 
as the counterparts of the turbulent MHD energy and cross helicity, Pk [Equa- 



tion (12a)] and Pw [Equation (13a)], respectively, with the reversed sign. This 
shows the productions of the turbulent MHD energy (u /2 + b' )/2 and cross he- 
licity (u ; • b') correspond to the sinks or drains of the large-scale MHD energy 
(U 2 + B 2 )/2 and cross helicity U • B. In this sense, the evolution of these turbulent 
quantities are subject to a kind of cascade process in turbulence. 



Appendix B. Turbulent cross- helicity expression based on the TSDIA 

In terms of the Elssasser variables 

</> = u + b, ift = u — b, 



(Bl) 



the magnetohydrodynamic (MHD) equations of the incompressible fluid can be 
written as 



50 
~dt 



+ (lf> ■ V)0 = -Vp M + ^V 2 0, 



(B2) 



V-0 = O, 



(B3) 



and the counterparts obtained by the interchange of <fi -<->• ip. Here, pm(= p + h 2 /2) 
is the MHD pressure. We have assumed the difference between the viscosity and 
magnetic diffusivities are small compared to the sum of them and neglect the 
former, putting v = r/. As for the extension of the Elssasser formalism to the 
compressible case, the reader is referred to [56] (and also [29] in the context of 
turbulence). 

In the two-scale direct-interaction approximation (TSDIA) formalism, we intro- 
duce two-scale variables: the fast and slow variables as 



£(=x), X(=<5x); r(=t), T(=5t). 
As this result, the spatial and time derivatives are written as 





V = V* + <5V X ; 



d . d 



(B4) 



(B5) 



This shows the property of derivative expansion with respect to the slow variables 



with 5 being the scale parameter. With the aid of Equation (B4), we separate the 
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fast variation of fluctuation from the slow variation of mean as 



2g 



f = F(K;T) + f'^,X;r,T). 



(B6) 



The Elssasser variables are divided into means and fluctuations as 



$ + cf)\ tp = * + ?//. 



(B7) 



We substitute Equation (B7) into Equations (B2) and (B3) in a rotating system 



with the angular velocity of flp. With Equations (B4)-(B6), the governing equation 



for the lowest-order fields in 5 is written in the wavenumber space as 



dr 



+ ^Vo a (k;r) 



iZ abc (k) J J «5(k - p - q)dpdq ^(P; r)C(q; r) 
-i(k • B)0' o a (k; r) - e^Z^(k) (0' o c (k; r)^ c (k; r)) 



(B8) 



and the t/> counterpart. Here, B is the mean magnetic field, D a (k)(= 5 —k a k /k 2 ) 
the projection operator, and Z a (k) = £; a D (k). Hereafter we suppress the depen- 
dence on the slow variables, X and T, except when necessary. In Equation (B8), 



the MHD-pressure term has been eliminated by using the Poisson equation for it. 
The equation for the first-order field is given as 
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-iZ abc (k) 1 1 5(k - p - q)dpdq $?(p; r)^(q; r) 

ya (k;r) 



-^ c (k)^(k;r)S-?4+^ L, " i) 



dX b DTi 



dxf 



+ie d6c O^ J D a(i (k) ' 



fe 2 " 



k" 



dX{ 





6 e (k;r) + V e (k;r)) 



+i6 rf ^^D-(k)^ (0 o c (k;r) + Vo C (k;r)) 

-i(k • B)^(k; r) - 6 d6c ^^(k) (^(k; r) + ^(k; r)) . 



(B9) 



Here, cj>' s is the first-order field satisfying the solenoidal condition [See Equa- 
tion (B13) later]. The ip counterpart can be obtained by exchanging 



4> o V, 3> «-► *, <f>' <r+ t/)' , Bo -B 



(BIO) 



in Equation ( |B9[). 
In Equation (B8) we consider the turbulence fields that is free from factors 



generating the anisotoropy of fluctuations such as the rotation and the magnetic 
field, and denote them the basic fields as cf)^ and ?/> B . They obey equations similar 
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to the homogeneous isotropic turbulence as 



0t 



+ ^V B a (k;r) 



-iz abc (k) J J 6(k - p - q)d P ciq $?(p; ^)0 B c (q; r) = 0. (Bll 

Corresponding to this equation, we introduce the Green's function 
dG'S b (k;r,r') 



„2rilab 



Or 



+ ^GT(k;r,r) 



■idbi 



tab i 



iZ acd (k) 8(k-p- q)dpdq V c (p; r)G7(q; r, r') = 5 ab 5(r - r'). (B12 



Using this Green's function, we formally solve the first-order field as 



^ 1 (k;r) = 0' sl (k;r)- i Ka<PB 



(k;r) 



with 



<Asi(k;r) 






lad/ 



k 2 dXf 



i lb I 



(B13) 



D cd (k) / driGn k ;^n)^(k;ri 



— oo 



d n GT(k;T, Tl 



D^(k;ri) 



DTi 



{(km) 



+B 6 / T (ir 1 G; flC (k;r,T 1 ) 

J— CO 

+ie dbc Q b F ^D d f(k)J T Jr 1 G f f(k; r, n)~ [^ B e (k; n ) + ^ e (k; n ) 



fc d 







+ie ^c^_ jD ae (k) / ^ G /a b(k;T)Ti) _^_ [^ (k; n) + ^ (k; n) ] . (B 1 4 ) 



The turbulent cross helicity is defined by 



(B15) 



Multiplying Equation (B13) with Equation (B14) by <//(k;T) and averaging, we 
integrate it over k. Then we obtain the expression for the turbulent cross helicity 

as 

W = 2I {Q ub } 

Iq "j Gs, — {Qub + Qbu) \ + h \ Ga, -jtj. (Quu + Qbb) 

+ ^ (O + 2ft F ) • [/_! {G s ,VH bu } + /_! {G A , VH uu }} 

-B • [/o {G s , V (Q UM + Q 6fe )} - /o {G A , V (Q Mfe + Q bu )}] , (B16) 



where use has been made of the abbreviated form of integrals [Equation (38)]. Here 
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we have introduced the statistically isotropic Green's function 

'Gt(k;T,T'))=6 ab G#(k;T,T f ) (B17) 



with $ = {(/>,ip), and defined the mirrorsymmetric and anti-mirrorsymmetric parts 
of Gff, and G$ , Gs and Ga , by 

G s (k; t, t') = (1/2) [G^k; r, t') + G^(k; r, r')] , (B18a) 

G A (k; r, t') = (1/2) [G^k; r, r') - G^k; r, r')] . (B18b) 

Appendix C. Suggestions from higher-order terms 

We consider situations where the time scales are mirrorsymmetric. In such a case, 
the C?A-related terms may be dropped. If we retain only the Gs-related terms in 



Equation (B16), W is expressed as 



W = 2I {Q ub } -iJgs,^ (Q ub + Q bu ) 

+ ^ (O + 2Q F ) ■ /_! {G s , VH bu } - B ■ J {G s , V (Q uu + Q bb )} . (CI) 

Here, Q uu , Q bb , and H bu are the spectral functions of the turbulent quantities at 
the lowest order, which are defined in terms of the basic or lowest-order fields u' B 
and bg as 

l - <u' B 2 + b' B 2 > = 2 f [Q uu (k; r, r) + Q bb (k; r, r)]dk, (C2) 



<b' B • u' B ) = f H bu (k; r, r)dk. (G3) 



In a manner similar to Section 4.2 we estimate the higher-order or the third and 
fourth terms in Equation (CI). We assume the spectral functions of the turbulent 
MHD energy and torsional correlation are written as 

Q uu (k,x;T,T',t) + Q bb (k,x;T,T,t) = a K (k,x;t)ex.p [-uj K (k,x; t)\r - r'\] , (C4) 

H bu (k,x;T,r',t) = a H (k,x]t)exp [-cj H (k,x.; t)\r - r'\] . (C5) 

Here, ok and an are the power spectra of the turbulent MHD energy and torsional 
correlation, respectively, and ujk and ujh reflect the time scales of fluctuations. (In 
this AppendixO the subscript H denotes the cross torsional correlation H bu not the 



residual helicity.) Note that in Equation (C4) we assumed the time-scale difference 



between the velocity and magnetic fluctuations is negligible and denote both the 
frequencies as ujk- The spectrum and frequency (reciprocal of time scale) of the 
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turbulent MHD energy can be expressed in the form 



(C6) 



u K (k, x; t) = uj K o£ 1/3 k 2/3 



' w 



(C7) 



(ctko, ^>ko are numerical factors). This corresponds to the Kolmogorov scaling of 
the turbulent MHD energy. 

We also assume the spectrum of the torsional correlation obeys a power law, and 
put 



a H (k,x.;t) = a m e 1/3 (x; t)e H (x; t)k~ 



11/3 



(C8) 



u H (k, x; t) = UHQe^k 1 ! 3 = r 



-i 
H ■ 



(C9) 



This may be a strong assumption, but as far as the spectral indices are concerned, 
they do not affect the final form of the dissipation-rate equation but only change the 
model constants. Alternatively, we may consider that a finite torsional correlation 
arises from the combination of the helicity and cross-helicity correlations such as 



a H {k, x; t) = a m e 4/3 (x; t)e H u{^\ *)ew(x; *) 



-11/3 



(CIO) 



(ehu- the helicity dissipation rate). In any case, we need further information on the 
torsional correlation or helicity spectrum, which so far we can not fully utilized. 



So, in this work, we assume the simple forms of Equations (C8) and (C9), and see 
the consequence. 



higher-order contributions in Equation (CI). Then we have 



We substitute Equations (C6) and (C8) with Equations (C7) and (C9) into the 



-B-I {G s ,V(Quu + Qbb)} 



B 



dk 



d,TiG s (k, x; r, n, t) V [Q uu (k, x; r, n, t) + Q b b(k, x; r, n, t)] 



-B a / dk 



1 3<Tx(^,x;t) crfc(k, x;t) dujK(k,x;t) 



US + U K 



dx a 



(us + uk) 2 dx a 



(Cll) 



(n + 2n F )-i_ 1 {G Sl VH bu } 



^(n a + 2n%) fk- 2 dkf 1 dT 1 G s (k, X ;T,T 1 ,t)VH bu (k,x-r,T 1 ,t) 



-(n a + 2n^) I k- 2 dk 



1 daH(k,x;t) o"i^(A;,x;t) duH(k,x;t) 
us + uh dx a (us + uh) 2 dx a 



(C12) 
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From Equations (C7) and (C9), we construct synthesized time scales as 



= 1 = oos(k,x;t) +u} K (k,x;t) 

7"SK t s T K 



(C13) 



111 /, ^ 

= 1 = ws(«, x; t) + uj H (k, x; t) 

TSH T S TH 



<^so 



1 urn / gg/fe \ 1/3 



W S0 



1/3^2/3 



^sh£ 



l/3 fc 2/3_ 



(C14) 



Namely, for the simplicity of argument, we assume that the time scales of turbulence 
are determined primarily by the energy transfer rate e. If we substitute these 



expressions (C13) and (C14) for ws + ojk and ujq, + ujh into Equations (Cll) and 



(C12), respectively, we have 



-B-I {G s ,V(Quu + Qbb)} 

= ^Tn—c 1/3 ^ 3 B a 
(2vr) 1 /3 Wsk C 



2 1uko\ 1 dg 

3 2 w sk / e <9x a 



11 + wko\ l_d£c 
3 ! w sk / ^c dx a 



(C15) 



-(fi + 2fi F )-/_i{G s ,Vif fe 4 

15(2vr) 7 / 3 u sh c v w \ edx a 

r 4 / m 1 ^g^ r j / ,,1 <9^C ] 

+ [3 - A h (uso,uho)]—-q-z + [11 + A H (uso,u m )] — — ^ 



(C16) 



with 



Ah(uso,uho) — 



30 u; g0 / gglcA 
33(2vr) 1 /3 Wsh V e ^ 



1/3 



(C17) 



Using an algebraic relation between £q and K: 



K = 3(2vr) 1 /V i , g 2/3 4 /3 



(C18a) 



or 



£ c = ^\2^)- 1 / 2 a K fe- 1 K^ 2 



(C18b) 



we express Equations (C15) and (C16) partly in terms of K instead of £q. Then 
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we have 



-B-I {G S ,V(Q UU + Q bb )} 

(2vr) 1 /3 Usk c 



3 I 2 uj &k ) e dx a 



1 /ll u K0 \ ^dK_ 
6 V 3 uj sk J K dx a 



(C19) 



(« + 2n F )-/-i{Gs,V# 6 u} 



1 



OffO 



15(2tt)V3 0J sh 



mf(^\[z c {w + 2^)] 



•[12 + ^ ff (u;so, Wifo)] £ ^ + [3 - Aj(wso,W£ro)] — ^§ 



+ 



33 3 

y + 2^( w so,who) 



1 3K 

Kdx" 



(C20) 



Equations (C19) and (C20) show that, in addition to the terms in Equation (62), 



the cross-helicity dissipation equation may contain the higher-order terms related 
to the mean magnetic field B and the mean vorticity fi (and the angular velocity 
ftp) such as 



K 3/2 



-(B-V)K, (B-V)e, K 1 / 2 (O • V) K, 

A £ 



(O-V)e. 



(C21) 



Model constants for these terms can be estimated by time scales of turbulence, uj^, 
u)Ki an d ^H- At the same time, inhomogeneity of the cross-torsional correlation 
(b ; • u/) , which is related to the coupling of the helicity and cross helicity, may play 
a certain role in the cross-helicity dissipation. 



